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SUMMARY 


The  thermomechanical  behavior  of  unidirectionally  fiber 
reinforced  laminae  in  plane  stress  has  been  described  by  a  set 
of  nonlinear  thermomechanical  stress-strain  relations.  A 
general  procedure  for  analysis  of  in-plane  loaded  symmetric 
laminates,  composed  of  such  laminae,  has  been  established  for 
simultaneous  loading  and  temperature  inputs.  A  ±45®  laminate 
subjected  to  heating  and  uniaxial  load  has  been  analyzed  in 
detail. 

Analyses  of  Graphite/Polyimide  laminates  have  shown  that 
the  thermomechanical  strains  cannot  be  separated  into  mechani¬ 
cal  strain  and  free  thermal  expansion  strain.  Consequently 
the  free  thermal  expansion  coefficient  is  generally  not  useful 
for  computation  of  thermal  strains  in  a  loaded  nonlinear 
laminate.  Moreover,  a  laminate  such  as  the  ±45®  which  behaves 
isotropically  for  free  thermal  expansion  behaves  anisotrop- 
ically  when  loaded  and  heated. 

Elastic  properties  and  thermal  expansion  coefficients  of 
unidirectional  Graphite/Polyimide  specimens  were  measured  as  a 
function  of  temperature  to  provide  inputs  for  the  analysis. 

The  ±45®  symmetric  Graphite/Polyimide  laminates  were  tested  to 
obtain  free  thermal  expansion  coefficients  and  thermal  expan¬ 
sion  coefficients  under  various  uniaxial  loads. 

The  experimental  results  demonstrated  the  effects  pre¬ 
dicted  by  the  analysis,  namely  dependence  of  thermal  expansion 
coefficients  on  load,  and  anisotropy  of  thermal  expansion 
under  load.  Numerical  agreement  of  analytical  and  experimen¬ 
tal  results  was  fair.  The  significance  of  time  dependence  on 
thermal  expansion  was  demonstrated  by  comparison  of  measured 
laminate  free  expansion  coefficients  with  and  without  15-day 
delay  at  intermediate  temperature. 


1.  INTRODUCTION 


Determination  of  the  thermoelastic  properties  of  fiber 
composite  laminates  on  the  basis  of  the  thermoelastic  proper¬ 
ties  of  the  constituting  laminae  and  their  stacking  sequence 
is  a  routine  procedure  used  by  composite  materials  analysts. 
Even  when  the  thermoelastic  properties  of  the  laminae  are 
temperature  dependent,  the  determination  of  laminate  proper¬ 
ties  is  of  conventional  nature,  provided  the  laminae  are 
elastic  at  all  temperatures. 

In  recent  years,  nonlinear  mechanical  properties  of  lam¬ 
inates  have  been  successfully  incorporated  into  laminate 
analysis,  primarily  in  the  isothermal  case.  In  the  case  of 
Boron/epoxy  and  Graphite/epoxy  layers,  room- temperature  non¬ 
linearity  is  exhibited  primarily  in  shear  and  transverse  ten¬ 
sion  and  compression.  Room- temperature  laminate  stiffness 
properties  are  affected  by  this  nonlinearity,  but  not  drasti¬ 
cally.  The  most  significant  nonlinearities  are  exhibited  by 
the  ±45°  laminate  under  uniaxial  loading. 

The  situation  appears  to  be  significantly  different  in 
the  case  of  laminates  subjected  to  large  temperature  changes. 
This  is  due  to  two  phenomena.  The  first  is  that  nonlinear 
behavior  of  a  unidirectional  lamina  results  from  matrix  non¬ 
linearity  which  increases  significantly  with  temperature. 

The  second  is  that  large  thermal  expansion  coefficients  exist 
in  directions  normal  to  the  fibers  of  a  unidirectional  lam¬ 
inae.  Therefore,  as  a  result  of  temperature  changes,  large 
transverse  stresses  or  strains  will  be  induced  even  for  a 
fiber-dominated  laminate.  It  is  to  be  expected  that  non¬ 
linearity  in  the  fiber  direction  can  be  neglected  since  the 
fibers  remain  stiff  and  elastic.  However,  nonlinear  effects 
of  the  stress-strain  relation  in  a  direction  transverse  to 
the  fibers,  and  in  shear,  both  increase  with  temperature  and 
have  to  be  considered.  Moreover,  these  transverse  and  shear 
effects  are  not  superposable,  but  are  interactive;  i.e.  the 
nonlinear  shear  strain  depends  also  on  the  transverse  strain 


2 


and  vice  versa.  These  nonlinearities  niay  have  a  very  signifi-* 
cant  effect  upon  thermal  expansion  coefficients  even  when  the 
effect  upon  laminate  stiffnesses  is  small. 

In  addition,  epoxies  and  polyimides  exhibit  time-dependent 
effects  whose  significance  also  increases  with  temperature. 
Therefore,  internal  stresses  and  strains  change  even  when  tem¬ 
peratures  and  loads  are  held  constant.  In  the  first  approxi¬ 
mation,  such  behavior  can  be  modeled  as  viscoelastic. 

The  above  effects  can  be  expected  to  occur  to  some  extent 
in  laminates  designed  for  moderately  elevated  temperature  ap¬ 
plications.  The  possible  consequences  of  these  effects  in¬ 
clude  the  following: 

1)  laminate  thermal  expansion  coefficients  and  internal 
stresses  as  determined  by  conventional  thermoelastic 
analysis  will  be  incorrect; 

2)  laminate  thermal  expansion  coefficients  and  internal 
stresses  will  be  time  dependent;  and 

3)  the  deformations  and  stresses  of  a  laminate  subjected 
to  loads  and  temperature  changes  will  not  be  correct¬ 
ly  obtained  by  adding  stresses  and  strains  due  to 
loads  alone  to  those  obtained  due  to  temperature 
changes  alone. 

The  temperature— dependence  problem  has  been  reported, 
e.g.  in  references  [1]  and  [2].  The  general  approach  is  that 
one  takes  as  a  starting  point  the  stress-free  temperature. 

(It  is  important  to  emphasize  that  depending  upon  the  fabrica¬ 
tion  process,  this  may  actually  mean  a  temperature  different 
from  room  temperature  at  which  all  stresses  in  the  laminate 
vanish.)  One  may  then  imagine  an  experiment  in  which  free 
thermal  expansion  of  each  layer  takes  place  from  the  stress- 
free  temperature  to  the  final  temperature  or  temperature  dis¬ 
tribution.  At  the  final  temperature,  one  may  then  apply 
stresses  to  satisfy  equilibrium  and  compatibility  and  reas¬ 
semble  the  laminate.  Hence,  it  follows  that  the  final  stress 
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distribution  is  obtained  by  using  secant  AL/L  values  from  the 
stress-free  temperature  to  the  final  temperature  and  by  using 
elastic  constants  at  the  final  temperature.  The  incremental 
theory  has  been  shown  in  reference  [2]  to  be  obtainable  from 
the  total  strain  theory  outlined  above. 

Time-dependent  properties  have  been  treated  by  a  number 
of  authors,  and  methods  of  varying  complexity  exist  in  the 
literature  (e.g.  refs.  [3]  to  [6]).  Similarly,  for  stress- 
dependent  phenomena  in  isothermal  laminate  analysis,  there  is 
a  significant  body  of  literature  (e.g.  refs.  [7]  to  [11]). 

For  both  of  these  problems,  the  various  existing  analytical 
methods  have  various  limitations. 

In  the  present  work  there  is  established  a  general  pro¬ 
cedure  to  compute  the  thermomechanical  properties  and  internal 
stresses  of  laminates  composed  of  nonlinear  temperature  depen¬ 
dent  laminae.  The  methods  used  are  similar  to  those  employed 
in  reference  [7]  for  isothermal  nonlinear  laminates. 

In  addition  to  the  analytical  investigation  an  experimen¬ 
tal  program  has  been  carried  out  whose  purpose  is  to  measure 
thermal  expansion  coefficients  of  free  and  loaded  Graphite/ 
Polyimide  laminae  in  order  to  examine  the  validity  of  the 
analytical  results. 


4 


2 .  THEORY 

2.1  Lamina  Thermomechanical  Stress-Strain  Relations 

The  approach  to  be  adopted  is  a  generalization  of  the 
method  of  reference  [7]  for  isothermal  nonlinear  stress-strain 
relations  for  plane  stress. 

The  lamina  state  of  stress  is  ^22'  '^12'  ^1 

in  the  fiber  direction  and  X2  is  in  the  transverse  direction 
(fig.  1)  . 

The  strains  in  the  lamina  plane  are  given  by: 


where  <|)  is  the  temperature. 

The  compliances  in  (2.1.2)  can  also  be  written  in  the 

form: 


5 


22  E^(<))) 


66  4G^((J)) 


(2.1.3) 


where  E  ,  v  ,  E  ,  G  are  temperature  dependent  axial  Young's 
modulus,  axial  Poisson's  ratio,  transverse  Young's  modulus, 
and  axial  shear  modulus,  respectively. 

It  is  assumed  that  because  of  large  fiber  stiffness,  the 

It 

nonlinear  strain  vanishes  and  that  has  no  effect  on 

II  i- 1  11  XI 

the  strains  ^12*  Consequently,  the  nonlinear  strains 

can  be  written  as: 


"ii  =  ° 


^22  “  ^22  ^*^22'  '^12'  ^22 


(2.1.4) 


^12  =  2Sgg  (0^2,  ^12 

M  11 

where  S„„  and  are  nonlinear  temperature  dependent  (secant) 

ZZ  o  b 

compliances • 

The  one-dimensional  isothermal  responses  to  0^2 
^12'  them  acting  separately,  are  modeled  in  the  famil¬ 

iar  Ramberg-Osgood  form: 


£ 

22  22 


E  +  £ 

12  ^  ^12 


^22  \ 


(2.1.5) 


where  and  are  temperature  dependent  stress  parameters 
and  m  and  n  are  exponents  which  may  also  be  temperature  de¬ 
pendent.  These  four  quantities  are  to  be  regarded  as  curve 
fitting  parameters  to  best  fit  the  experimental  data. 


6 


Employing  the  same  reasoning  as  in  [7] ,  it  follows  that 
in  the  case  of  combined  stresses  ,  the  nonlinear  parts  of 
the  strains  are  given  by: 


Finally,  the  free  thermal  strains  in  (2.1.1)  are  written 

in  the  form: 


(<|))<|) 


^22  “t 


(2.1.7) 


where  a  and  a  are  axial  and  transverse  temperature-dependent 

J\  J. 

(secant)  thermal  expansion  coefficients  assumed  to  be  stress 
independent . 

Summarizing  equations  (2.1,1)  through  (2*1,7),  we  have 
the  following  nonlinear  temperature  dependent  stress-strain 
relations: 

Hi  =  i^“ll  -  i^"22  *  “a 


_  1 
^22  “■■E^(<j))  *^11  ^  E^(4)) 


ay(<|.) 


Ty(<}.) 
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2,2  Analysis  of  Nonlinear  Temperature  Dependent  Laminates 

Consider  a  balanced  symmetric  laminate  loaded  on  its 
edges  by  uniformly  distributed  stresses  and  subjected  to  the 
temperature  rise  4). 

The  following  problems  are  to  be  considered:  (i)  find 
the  stresses  in  the  laminae;  (ii)  find  the  effective  thermo¬ 
mechanical  stress-strain  relations  of  the  laminate. 

The  laminate  is  referred  to  a  fixed  system  of  coordin¬ 
ates  ,  x^  (fig*  1) •  Any  lamina,  kth  say,  is  referred  to 

.  (k) 

its  material  system  of  coordinates,  x^  in  fiber  direction 

(k)  ^ 

and  x^  transverse  to  fiber  direction  (fig.  2) .  The  fiber 

orientation  of  the  kth  lamina  is  defined  by  the  angle: 


(2.2.1) 


The  stress-strain  relation  of  the  lamina  in  its  material 


system  is 


S  a  +  a 


2.2.2) 


where  underbars  denote  tensors.  Here  is  the  temperature  and 
stress  dependent  compliance  tensor  and  a  the  temperature  de¬ 
pendent  thermal  expansion  tensor  as  defined  by  (2.1.8). 

It  may  be  easily  shown  by  trivial  considerations  of  heat 
conduction  theory  that  the  temperature  throughout  all  laminae 
is  the  same  constant  (}) .  Also  in  accordance  with  usual  lamin¬ 
ate  theory:  (i)  all  laminae  are  in  states  of  plane  stress; 

(ii)  the  strains  in  all  laminae  as  referred  to  the  laminate 
fixed  coordinate  system  are  the  same  uniform  . 

These  considerations  make  it  possible  to  construct  a  set 

of  nonlinear  equations  for  computation  of  all  laminae  stresses. 

(k) 

First,  define  the  transformation  tensor  T  by: 


(k)^  ^  T(k)^(k) 


2.2.3) 
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where  a is  the  lamina  stress  referred  to  the  material  sys¬ 
tem  and  is  the  same  stress  referred  to  the  laminate 

system. 

(k) 

The  matrix  of  T  is  given  by: 


'Cos^Gj^  Sin^Gj,^  -2CosGj^Sin6j^' 

r>p(k)j  _  sin^Gj^  Cos^Gj^  2CosG  SinGj^ 

^  2.2 

Sin0j^Cos6j^  “Sin9j^Cos9j^  Cos  9^-Sin  9j^ 

(2.2.4) 


Similarly, 

(k),  .  jk  ^(k) 


(2.2.5) 


It  follows  from  equilibrium  that 

^  =  °ah  (2.2.6) 

1^1 

where : 

2K  -  number  of  laminae; 

2h  -  thickness  of  laminate; 
tj^  -  thickness  of  kth  lamina; 

o  o  o  o 

and  a  has  components  the  stresses  *^22' 

aj^2  applied  to  the  laminate  edges. 

By  hypothesis  (b) 

(1),  _  (2),  _  _  (k)e  _  (K)^  (2.2.7) 

Now  introduce  (2.2 .3)  into  (2.2.6),  (2.2.2)  into  (2.2.5) 
and  subsequently  into  (2.2.7).  The  results  are: 
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II  ^ 


—  k  — 


(2.2.8 


(1).  ,(1)  ^  md) 


(2)^  ^(2)  ^  „(2) 


o'  '  +  T'-'a 


o‘’'>  +  T'^’a-f  . 


T  (  k)  s  0  _j_  jp  (k) 


(2.2.9) 


O 

since  a  has  three  components,  there  are  altogether  3K  un¬ 
known  laminae  stresses.  To  determine  these  we  have  the  three 
equations  (2.2.8)  and  the  3(K-1)  equations  (2.2.9),  in  all  3K 
nonlinear  equations . 

Once  the  laminae  stresses  are  known,  the  effective 
thermomechanical  stress-strain  relations  of  the  laminate  can 
be  determined  in  the  following  fashion.  Write  (2.2.2)  in  the 
laminate  coordinate  system: 


e  =  S  a  +  a 


2.2.10) 


All  tensors  in  (2.2.10)  are  obtained  from  their  counterparts 
in  (2.2.2)  by  the  transformation  (2.2.3). 

Define : 


(k)c  =  (k)g-l 


2.2.11 


where  the  inversion  has  to  be  carried  out  for  components 

computed  as  functions  of  the  lamina  stresses.  Since  also  by 
(2.2.7)  ^  is  the  laminate  uniform  strain  £,  (2.2.10)  can 

be  written  as: 


(k)^  =  _  (k)^(k)^ 


(2.2.12 
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Substituting  (2.2.12)  into  (2.2.8)  gives; 

0*0* 

a  =  C  e  +  r  <{) 


(2.2.13) 


K 

r*  =  -  f  t,/h  (b)  (2.2.14) 

k=l  -  -  ^ 

Equation  (2.2.13)  formally  resembles  an  effective  thermo¬ 
elastic  stress-strain  relation.  It  is  nonlinear  since  C  and 
* 

£  depend  on  the  laminae  stresses. 

The  inversion  of  (2.2.13)  is; 


0  *  o  * 

£  =  S_  £  +  a  (|i 

where 

*  *  -1  *  0 

S  =  (C  )  =  S  (a  ,4)) 

(a) 

*  *  *  *  O 

a  =  -S  r  =  a  (a  ,<()) 

(b) 

(2,2.15) 


(2.2.16) 


It  is  seen  that  (2.2.16a)  are  the  nonlinear  compliances  of  the 
laminate  while  (2.2.16b)  are  its  nonlinear  stress  dependent 
thermal  expansion  coefficients. 

The  usual  definition  of  thermal  expansion  coefficients  is 
based  on  heating  of  an  unloaded  body.  If  this  is  done  in  the 
present  case,  equations  (2.2.6)  -  (2.2.8)  become; 


k=l 


K 

Z  -0 

k=l  ^ 


(2.2,17) 
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The  last  of  (2.2.17)  together  with  (2.2.9)  now  determine 
the  laminae  stresses. 

Repeating  the  same  procedure  as  previously  used  it 
follows  that: 

£  =  a  cj)  (2.2.18) 

where  a  is  defined  by  (2.2.16b),  (2.2.14b)  and  (2.2.11). 

Note,  however,  that  since  the  laminae  stresses  in  the  present 
free  body  heating  problem  will  be  different  from  the  previous 
ones,  the  thermal  expansion  coefficients  defined  by  (2.2.18) 
are  in  general  different  from  the  ones  obtained  in  the  case 
of  laminate  loading  and  heating.  Indeed,  the  thermal  expan¬ 
sion  coefficient  in  (2.2.18)  can  be  written  as  a  special  case 
of  (2.2.16b);  i . e . : 


a  =  ot  (0 ,  (j))  (2.2.19) 

("a=0) 

Consequently,  the  classical  concept  of  thermal  expansion 
coefficient  is  not  useful  in  the  case  of  a  nonlinear  laminate 
since  the  strains  due  to  loading  and  heating  can  no  longer  be 
found  by  super-position  of  mechanical  and  thermal  strains. 
Instead,  the  laminate  strains  must  be  found  on  the  basis  of 
the  interactive  effect  of  loading  and  heating. 


2.3  Analysis  of  ±45°  Laminate 

The  general  analysis  of  paragraph  2.2  will  now  be  applied 
to  the  case  of  a  balanced  ±45°  laminate  which  is  heated  and 
subjected  to  uniaxial  stress  (fig.  3).  From  symmetry  con¬ 
siderations: 


a 


1 

11 


.3.1) 
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The  stress  transformation  (2. 2. 3-4)  with  6  =  45°  yields: 


1  _  1  1  ^ 

2  2  ^22  ■  12 


11.11  1 

2  ‘^ll  2  ^22  '^12 


I 

2  11  2  22 


(2.3.2) 


The  three  equations  of  equilibrium  (2. 2 ,6)  assume  the 


form: 


12  ° 

°ii  '^ii  ^  ^  '^ii 


^^22  ^^22  ° 


1^2  . 
012  +  ^12  =  0 


(2.3.3) 


In  view  of  (2.3.1)  and  (2.3.2),  equations  (2.2.3)  become: 


”u  * 


“ll  *  <’22  ^  2”i2  =  “ 


(2.3.4) 


the  last  of  equation  (2.3.3)  is  identically  satisfied 
Adding  and  subtracting  equations  (2.3.4)  gives: 


11° 

*^11  '*■  ^22  *^11 


1  ° 

2  ^11 


(2.3.5) 


thus  determining  the  lamina  shear  stress 
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The  lamina  strains  are  given  by 


"n  +  +  ^12^22  ^ 


^22  ^12^11  ^22*^22  ■'■  “2'*' 


where  according  to  (2.1.8): 


11 


12 


EA(ct>) 

Ea(^) 


22  E^((J)) 


G  Ga(<J>) 


1  + 


1  + 


(aj^))  ^  (t 


12 


(cf.) 


i-Z22\ 

la  (<t>)J 


U,.  {<!))/ 


“1  "" 


“2  “t^***^ 


From  symmetry 


^^12  ^12  ° 


Therefore : 


^11  ^22 


(2.3.6) 


¥ 


(2.3.7) 


(2.3.8) 


(2.3.9) 
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Introducing  (2 . 3 . 6)  into  (2,3.9)  gives: 


=  ir=12>  <=12-=22>  “22 


2.3.10 


Equation  (2.3.10)  contains  the  stresses  0^2  ^ 12  non¬ 

linear  form.  Together  with  (2.3.5)  it  determines  the  stresses 
and  022*  Eliminating  between  these  two  equations  and 

using  (2.3.5b)  yields  the  result: 


1  \  ^ 

^^22'  ^  ^  ^  1,^  ^ 


1  2 


2-,  m-l- 


<1  +  “a>  “11  *  *  0 


2.3.11) 


which  determines  ^^22* 

The  laminate  normal  strains  are  given  in  terms  of  tensor 
transformation  and  use  of  (2.3,9)  by; 


.  - 

^11  11  11  12 


1  -1^1 
^22  “  ^22  ^11  ^12 


(2.3.12) 


From  (2.3.7),  (2.3.6),  (2.3.5a)  and  (2.3.12 
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O  1 

^22  ^  ®22 

o  -  _  1  2  o 


(  In  ' 

I  o..  ^ 


“a^ 


2.3.13 


Next  some  special  cases  are  considered:  Let  the  laminate 
first  be  heated  without  loading.  Then  (2.3.11)  becomes: 


1  m-1 


1  ^  2“a  i; 


Jr  ®a'W'»  °  “ 

i-  V  I 


2.3.14) 


The  laminate  strains  are  then  from  (2.3.13): 


£11(0, 4*)  =  £22^®''^^  “ 


1+^A  1 

ST  ”22  V 

A 


(2.3.15) 


where  in  (2.3.15)  is  defined  by  (2.3.14).  Consequently, 
the  free  thermal  expansion  coefficient  of  the  laminate  is 
given  by: 


*  ^"^^A  1 

“f  ^  ^22  “a 


(2.3.16) 


Now  let  the  laminate  be  subjected  to  the  axial  stress 

o  ]_ 

with  ct)=0.  Then  0^2  defined  from  (2.3.11)  by: 


«  2.  m-1 

^ll)  ~ 
2t 


^  ^A^  '^11 


2.3.17) 


and  the  strains  are: 
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1°  0\  =  ^11  1 

22 


1  2 


2i  n-1 


"22  (  ^11,0^ 


^  1 


,  21  n-1 
^11  \  2 
2t_. 


(2.3.18 


By  conventional  definition  the  expansion  coefficients  are 
given  by: 


“l  (p  "ll^  '^11,'*’^  ~  "ll^*^!!,®^ 


*  1  ®  /  rr  ®  ® 

“2  ~  ^  ^22  “  "22^'^ll,®^ 


(2.3.19) 


It  is  seen  that  the  expansion  coefficients  as  defined  by 
(2.3.19)  are  load  dependent,  different  from  the  free  expan¬ 
sion  coefficients  (2.3.16)  and  in  fact  unequal .  Thus  the 
thermal  expansion  of  the  laminate  has  become  anisotropic  due 
to  load  application. 

Conventional  linear  thermoelastic  analysis  is  contained 
in  the  preceding  results  as  a  special  case  when: 


°y  '^22 


T  >>  a,- 

y  12 


which  implies  that  nonlinear  effects  are  negligible.  In  this 
case  we  have  for  free  thermal  expansion  from  (2.3.14)  and 
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2.3.16)  : 


1 

22  l+2v^+E^/E^ 


“f  ~  “a  1+2v^+E^/E^' 


(2.3.20) 


which  are  the  well-known  thermoelastic  results. 

It  also  follows  from  (2.3.11)  and  (2.3.13)  that  for 
heating  and  loading; 


eil(  C7^j_,(f)) 


11  * 


o  '^12®  * 

^22<  Oil.^)  =  -  — 5-  '^11  +  “f*^ 

^1 


(2.3.21) 


"ic  "fc  '  , 

where  E^  and  v^2  laminate  Young's  modulus  and  associ¬ 

ated  Poisson's  ratio  at  temperature  (|) . 
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3,  EXPERIMENTAL  PROCEDURE 


Strain  Gage  Techniques 

For  the  constant  temperature  mechanical  tests  the  quarter 
bridge  strain  gage  circuit  was  used  since  no  temperature  com¬ 
pensation  was  needed.  A  temperature  compensated  half  wheat- 
stone  bridge  circuit  was  used  for  measuring  the  material 
thermal  response  over  the  -53®C  to  260®C  (-65®F  to  +500®F) 
temperature  range  to  eliminate  the  nonlinear  effects  from  the 
changes  in  resistivity  of  the  gage. 

In  the  half  bridge  circuit  (fig.  4)  the  active  gage  was 
bonded  to  the  test  specimen  and  the  compensating  gage  bonded 
to  the  quartz  standard  with  a  known  coefficient  of  thermal 
expansion  (Dynasil  Corp.,  a  =  1.82ye/°C).  The  resulting 
strain  output  from  this  circuit  is  the  difference  between  the 
active  and  compensating  gage.  Knowing  that  the  relationship 
between  the  apparent  strain  and  the  thermal  coefficients  of 
the  gage  and  substrate  is: 


a  =  a  +  e  /Acf) 
s  g  app 

where ; 


a  =  coefficient  of  thermal  expansion  of  substrate; 
s 

=  coefficient  of  thermal  expansion  of  grid  (gage) ; 

e  =  apparent  strain; 

app 

Y  =  coefficient  of  resistivity  of  gage; 

Sg  =  gage  factor. 

For  a  given  temperature  excursion 

a  =  a  +  e  /A<()  -  y/S  ;  (3.1) 

p  g  P  g 

a  =  a  +  e  /A(|)  -  y/S  ?  (3.2) 

q  g  q  ^  ^  g 


where  subscripts  p  and  q  refer  to  the  Graphite/Polyimide 
substrate  and  quartz,  respectively. 
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Subtracting  and  assuming  identical  gages  on  each  material: 


a 

P 


a  = 

q 


e 


where  the  resistivity  and  a  grid  are  eliminated.  The  appar¬ 
ent  strain  difference  (e  -  e  )  is  measured  in  the  half 

P  q 

bridge  circuit.  Therefore ,  the  thermal  coefficient  of  expan¬ 
sion  for  the  polyimide  is: 


a  =  (e  -  G  )/A<()  +  a  .  (3.3) 

P  P  q  q 


Where  (g^  -  Gg)/A(j)  is  easily  found  from  the  slope  of  the 
measured  strain  versus  temperature  curve  and  is  known, 
1.82yG/°C  (1 . 01|iG/°F)  .  For  example,  the  measured  slope  of  the 
unidirectional  lamina,  a,  was  26.66yG/°C  (14.81yG/°F)  in  the 
transverse  fiber  direction,  making  a  =  28.48yG/°C  (15.82yG/ 
“F). 


The  gages  used  in  all  of  the  strain  measurements  were 
fully  encapsulated  ”K  alloy"  with  operating  temperature  ranges 
from  “269°C  to  +288°C  (-452°F  to  +550®F) .  Two  types  of  gage 
were  used,  one  having  a  thermal  coefficient  of  expansion  of 
6.0]jg/°C  and  the  other  having  a  coefficient  of  23,4ijg/®C.  The 
coefficient  of  expansion  values  were  chosen  to  maximize  the 
apparent  strain  output  by  taking  advantage  of  thermal  mis¬ 
matches.  The  apparent  strain  output  is: 


G 

app 


a  )  A(|)  +  yA(()/S  . 
y  y 


when  is  small  the  output  can  be  increased  by  making  dif¬ 
ferent  from  a  . 

s 

Temperature  measurements  were  made  using  a  bondable  resis¬ 
tance  temperature  sensor.  These  gages  were  used  in  combination 
with  a  Linearized  Strain  Temperature  (LST)  network  which  con¬ 
verted  the  temperature  reading  to  a  calibrated  strain  reading 
{18yG/®C) .  These  gages  were  bonded  to  each  sample  to  give  a 
temperature  reading  of  sample  surface  temperature. 
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An  epoxy  based  high  temperature  adhesive  with  an  opera¬ 
ting  temperature  range  from  -269°C  to  +371°C  (-452°F  to  + 
700®F)  was  used  to  bond  the  gages  in  place.  The  adhesive  was 
cured  at  a  temperature  of  163°C  (325°F)  for  two  hours  and  no 
post  cure  was  used.  All  samples  were  cured  at  the  same  tem¬ 
perature  and  for  the  same  period  of  time. 

At  all  times  during  the  pretest  period  the  samples  were 
stored  in  sealed  plastic  bags  containing  a  dessicant  to  absorb 
all  entrapped  moisture.  After  the  completion  of  all  wet 
machining  and  NDT  procedure  all  samples  were  dried  for  24 
hours  at  121°C  (250°F)  then  again  stored  in  a  sealed  plastic 
bag  with  dessicant. 

Mechanical  Characterization  of  the  Lamina 

Tensile  tests  were  conducted  at  room  temperature  -54 
(-65°F)  and  260°C  (500°F)  to  determine  the  dependence  of 
mechanical  properties  on  temperature.  The  tests  were  con¬ 
ducted  in  an  environmental  chamber  which  has  temperature 
control  from  -73°C  to  260°C  (-100°F  to  500°F) .  Air  is  circu¬ 
lated  constantly  throughout  the  test  chamber  to  maintain 
reasonably  uniform  temperatures.  Cooling  was  accomplished  by 
liquid  nitrogen.  Each  sample  was  held  at  the  test  temperature 
for  30  minutes  to  insure  that  equilibrium  temperature  had  been 
reached. 

The  0®  tensile  coupons  measured  1.27  cm  x  15.24  cm 
(1/2"  X  6")  and  the  90°  coupons  measured  2.54  cm  x  15.24  cm 
(1"  X  6").  End  tabs  were  manufactured  from  3.78  mm  (1/8") 
polyimide  circuit  board  material  which  was  cut  into  3.8  cm 
(1.5")  strips  and  beveled.  The  end  tabs  were  bonded  to  the 
coupons  in  a  special  fixture  to  assure  proper  alignment.  Each 
coupon  was  then  instrumented  with  a  longitudinal  and  trans¬ 
verse  strain  gage  (the  06  gage  was  used) . 

Chamber  temperature  was  measured  by  a  temperature  sensor 
mounted  on  a  dummy  coupon  placed  in  close  proximity  to  the 
coupon  under  test.  Strain  data  was  taken  at  known  load 
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increments  to  generate  stress-strain  curves.  From  these 
curves  E  ,  E  ,  v  ,  v  could  be  determined.  Approximately  20 

A  J.  J.r\ 

data  points  were  taken  from  each  sample;  three  samples  were 
tested  for  each  property  and  test  temperature.  A  linear 
least  squares  fit  was  used  to  determine  the  above  properties. 

Thermal  Expansion 

Thermal  coefficients  of  expansion  were  measured  for  uni¬ 
directional  lamina  and  [±45° ]  laminates.  Load  and  relaxation 

s 

effects  upon  the  coefficient  of  expansion  were  also  tested  on 

the  [±45° ]  laminate - 
s 

The  5  cm  X  5  cm  (2"  x  2")  test  specimens  were  suspended 
in  the  center  of  the  test  section  of  the  environmental  chamber 
for  all  tests.  Both  transverse  and  longitudinal  strain  gages 
were  mounted  on  the  unidirectional  samples,  while  a  single 
gage  was  sufficient  for  the  [±45°]^  laminates.  All  instrumen¬ 
tation  was  clustered  near  the  sample  center  to  avoid  errors 
due  to  edge  effects. 

For  both  the  unidirectional  and  [±45° specimens,  the 
temperature  excursion  was  -53°C  to  260°C  to  -53°C  (-65°F  to 
500°F  to  -65°F) .  By  a  complete  temperature  cycle,  hysteresis 
effects  could  also  be  observed.  Strain  data  were  recorded  at 
-13.9°C  (25°F)  increments  throughout  the  excursion.  A  plot  of 
the  strain  versus  temperature  for  each  sample  allowed  deter¬ 
mination  of  the  -  e_)/A(j)  slope. 

P  Q 

To  determine  the  variation  in  longitudinal  and  transverse 
thermal  coefficients  with  temperature  of  the  unidirectional 
laminae,  the  curves  were  partitioned  into  temperature  inter¬ 
vals  and  a  linear  least  squares  fit  applied  to  data  in  each 
interval  (the  data  in  the  intervals  were  essentially  linear) . 
The  coefficients  of  expansion  were  determined  in  each  interval 
by  applying  equation  3.3.  For  example,  in  the  -53°C  to  38°C 
(-65°F  to  100°F)  temperature  range,  for  sample  All, 

(e  -  e  )/A(|)  was  24.84ye/°C  (13.80yG/°F)  which  yields  an 
P  q 

a  =  26.66y£/°C  ( 1 4 . 81y£/°F) .  Similarly  in  that  same 
P 
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temperature  range,  for  sample  All,  (e  -  e  )  /A(()  was 

H  SI 

0.70y£/°C  (0.39yE/'’F)  in  the  transverse  direction;  which  makes 
=  1.12yE/°C  {0.62yE/'’F)  . 

The  [±45°],.  test  data  were  collected  in  the  same  manner 
as  above.  The  analysis  of  the  [±45]  data  consisted  of  a 
linear  least  squares  fit  applied  to  a  linear  portion  of  the 
curve  in  the  lO^C  to  149°C  (50°F  to  300°F)  temperature  range 
only. 

Compression  samples  were  constructed  for  the  loaded  tests. 
Each  sample  was  made  1.9  cm  (3/4")  in  width  with  a  1.58  cm 
(5/8")  test  section.  End  tabs  were  bonded  on  with  a  high 
strength  epoxy  adhesive  and  tape  was  used  to  prevent  bonding 
in  the  area  adjacent  to  the  test  section  thus  increasing  the 
test  section  to  one  inch  in  length  (2.54  cm)  but  with  only 
5/8"  (1.58  cm)  supported.  The  IITRI  compression  fixture  was 

used  and  the  constant  loading  was  accomplished  by  use  of  a 
closed  loop  servo-hydraulic  test  machine.  Tests  were  run  at 
150  kg  (5.13  ksi  or  35.37  MPa),  300  kg  (10.26  ksi  or  70.74  MPa) 
and  400  kg  (13.33  ksi  or  91.77  MPa).  Gages  were  mounted  both 
transversely  and  longitudinally  with  respect  to  the  loading 
direction.  Strain  was  recorded  at  13.9°C  (25°F)  temperature 
intervals  over  the  full  -53°C  to  260°C  to  -53°C  (-65°F  to 
500°F  to  -65°F)  cycle.  The  results  were  analyzed  in  the  same 
manner  as  the  other  [±45]^  results. 

Two  types  of  thermal  cycles  were  run  to  assess  relaxation 
effects.  One  group  of  samples  was  taken  from  260°C  to  room 
temperature,  stored  for  15  days  at  room  temperature,  then 
cooled  to  -53°C.  The  other  samples  were  taken  from  260°C  to 
room  temperature,  delayed  15  days  and  then  heated  to  260°C. 
During  the  delay  period  the  samples  were  stored  in  dessicant 
except  for  the  4.5  day  period  when  they  were  stored  in  an 
oven  at  121°C.  The  first  five  days  of  the  delay  period, 
samples  were  stored  at  room  temperature,  then  the  4,5  day 
121°C  treatment  was  administered  followed  by  a  5.5  day  room 
temperature  storage  period.  Strain  data  were  recorded  and 
analyzed  as  previously  described. 
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Test  Procedure  for  Thermal  Expansion  Coefficients 
of  Loaded  [±45]^  Laminates 

The  instrumented  sample  was  placed  in  the  compression 
fixture  which  was  centered  in  the  environmental  chamber.  In 
the  no-load  condition,  the  strain  gages  were  zeroed  and  load 
was  applied  to  the  sample.  An  initial  strain  reading  was  re¬ 
corded  to  indicate  the  no-load  strain  component  and  then  the 
chamber  temperature  was  reduced  to  -54°C  (-65°F) .  The  sample 
was  heated  to  260°C  (500°F)  at  a  rate  of  approximately 
83®C/hr.  with  strain  data  recorded  at  13.9°C  (25°F)  tempera¬ 
ture  intervals.  Cooling  was  initiated  immediately  at  260°C 
(500°F)  and  strain  data  were  recorded  at  13.9°C  (25°F)  inter¬ 
vals  for  the  descending  temperature  excursion  to  -54°C 
(-65°F) .  Total  test  time  was  approximately  seven  hours  for 
each  test. 

Figures  5-8  show  plots  of  longitudinal  and  transverse 
strain  (plotted  is  which  has  thus  to  be  adjusted  by 

the  quartz  thermal  expansion)  for  ascending  and  descending 
temperature  excursions.  All  of  the  plots  show  hysteresis 
effects  which  are  most  significant  in  the  case  of  longitudinal 
expansion.  The  large  increases  in  compressive  longitudinal 
strain  from  149°C  (300®F)  to  260°C  (500°F)  indicate  a  consid¬ 
erable  loss  of  stiffness  of  the  laminate  in  this  temperature 
range,  which  can  probably  be  due  to  large  creep  at  those 
temperatures.  The  stiffness  is  regained  upon  cooling  down. 
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4.  RESULTS  AND  DISCUSSION 


Test  results  have  been  obtained  for  the  temperature  de¬ 
pendent  elastic  moduli  and  thermal  expansion  coefficients  for 
the  unidirectional  material  and  for  ±45°  laminates.  The  ±45° 
laminates  were  selected  as  a  configuration  for  which  compari¬ 
sons  could  be  made  between  analytical  predictions  and  experi¬ 
mental  results.  These  laminates  were  examined  to  assess  the 
effects  of  nonlinear  stress-strain  behavior^  time  dependent 
stress-strain  behavior,  and  of  externally  applied  loads  upon 
thermal  expansion  coefficients  of  a  laminate. 

Properties  of  Unidirectional  Composites 

As  shown  in  table  1,  the  elastic  properties  of  the  uni¬ 
directional  material  are  accurately  reproducible  and  their 
variation  with  temperature  is  not  very  significant.  Results 
for  the  two  thermal  expansion  coefficients,  as  shown  in  table 
2,  show  somewhat  erratic  behavior,  both  with  respect  to  vari¬ 
ation  with  temperature  and  choice  of  specimen.  The  scatter 
is  more  pronounced  for  a  than  for  a  .  The  reason  for  this 

X 

is  that  a  is  a  much  smaller  number,  of  the  order  2ye°C 
and  the  accuracy  of  thermal  expansion  coefficient  measure¬ 
ments  for  such  small  values  is  questionable.  This  is  an  in¬ 
trinsic  problem  for  all  thermal  expansion  measurements  with 
conventional  apparatus . 

Properties  of  ±45°  Laminates 

A  number  of  different  experiments  have  been  carried  out 
on  these  ±45°  laminates.  The  major  variables  examined  were 
the  temperature  range  over  which  the  expansion  was  measured 
and  the  presence  or  the  absence  of  load  during  the  thermal 
expansion  test.  Experimental  results  are  presented  in 
tables  3  and  4.  In  order  to  attempt  correlation  between  the 
theoretical  and  the  experimental  results,  computation  of 
thermal  expansion  coefficients  has  been  carried  out  on  the 
basis  of  the  theory  presented  earlier  utilizing  measured 
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lamina  properties.  The  results  for  these  calculations  are 
shown  in  table  5.  Here  calculated  average  thermal  expansions 
are  presented  for  different  starting  and  finishing  tempera¬ 
tures.  Calculations  of  the  free  thermal  expansion  coeffi¬ 
cient  depend  upon  the  Ramberg-Osgood  parameter,  for  the 

transverse  stress-strain  relation  of  the  unidirectional  mate¬ 
rial.  This  value  is  not  well  known,  and  hence  the  computa¬ 
tions  have  been  carried  out  for  various  assumed  yield  stress 
values.  Decrease  of  the  value  of  a  implies  the  increase  in 

y 

nonlinearity.  An  infinite  value  of  implies  elastic  be¬ 
havior.  It  is  seen  that  the  effect  of  nonlinearity  does  not 
become  pronounced  unless  has  a  value  below  69  MPa. 

Free  Thermal  Expansion 

Experimental  results  for  thermal  expansion  of  the  ±45° 
laminate  between  10°  and  149 °C  is  shown  in  table  3a.  Compar¬ 
ison  of  the  mean  values  from  the  three  samples  tested  with 
the  computed  thermal  expansion  coefficients  for  the  same  tem¬ 
perature  range  as  shown  in  table  5  yield  only  fair  agreement. 
This  comparison  between  experimental  and  analytical  results 
is  not  good  enough  to  assess  the  importance  of  nonlinearity 
for  the  free  expansion  coefficients.  It  should  be  noted  that 
when  the  lamina  behavior  is  elastic  (whether  it  be  linear  or 
nonlinear)  the  average  thermal  expansion  coefficients  between 
two  temperatures  is  the  same  for  both  heating  and  cooling 
over  that  temperature  range.  The  experimental  results  indi¬ 
cate  a  difference  between  measured  thermal  expansion  on  the 
ascending  and  descending  sections  of  the  measurements. 

One  underlying  difficulty  associated  with  the  calcula¬ 
tions  is  the  need  to  assess  the  temperature  at  which  the  lam¬ 
inate  is  stress  free,  for  expansion  coefficient  calculations 
are  always  performed  relative  to  the  stress  free  state.  It 
is  reasonable  to  assume  that  a  laminate  which  has  been  kept 
for  sufficient  time  at  room  temperature  is  stress  free  since 
residual  stresses  will  diminish  with  time  by  a  relaxation 
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process.  When  the  laminate  is  heated  from  this  state,  inter¬ 
nal  stresses  develop.  However,  at  elevated  temperatures  such 
stresses  may  become  very  small  after  a  relatively  short  time 
since  stress  relaxation  is  much  more  predominant  at  elevated 
temperatures . 

Table  3b  illustrates  this  effect.  Here,  during  the  cool¬ 
down  from  260‘’C  to  room  temperature  it  is  reasonable  to  assume 
that  this  high  initial  temperature  was  the  stress  free  tem¬ 
perature.  The  samples  were  then  held  for  15  days  at  room  tem¬ 
perature  providing  an  opportunity  for  stress  relaxation.  At 
this  point  the  same  temperature  range  was  traversed  by  heat¬ 
ing  and,  as  can  be  seen  in  table  3b,  the  resulting  thermal 
expansion  coefficient  was  significantly  higher  than  during 
the  cooling  temperature  cycle.  Examination  of  table  5,  for 
the  temperature  range  of  room  temperature  to  260 °C,  shows  that 
there  is  reasonable  agreement  between  analysis  and  experiment 
for  the  case  in  which  260‘’C  is  taken  as  the  stress  free  tem¬ 
perature.  Calculations  for  the  heating  cycle  could  not  be 
made  because  of  the  uncertainty  of  material  elastic  proper¬ 
ties  at  the  final  high  temperature.  Such  properties  are  re¬ 
quired  for  the  final  temperature  of  the  heating  experiment 
as  discussed  earlier.  For  this  reason  subsequent  experiments 
were  conducted  at  lower  temperatures. 

The  effects  of  delay  time  were  treated  at  temperatures 
below  room  temperature  as  shown  in  tables  3c  and  3d.  Laminate 
expansion  coefficients  were  measured  after  cooling  from  260°C 
in  the  range  RT  to  -54°C.  Measurements  were  made  both  in  the 
case  of  a  15  day  delay  at  room  temperature  and  for  no  delay 
at  room  temperature.  It  is  assumed  that  initially  the  lamin¬ 
ate  was  stress  free  at  260°C.  When  the  laminate  is  cooled 
without  delay  to  -54°C  internal  stresses  develop  and  grow 
with  decrease  in  temperature.  In  the  case  of  15  day  delay  at 
room  temperature  it  is  assumed  that  the  stresses  which  have 
formed  in  cooling  from  260°C  to  room  temperature  relax  to 
such  an  extent  that  they  can  be  neglected.  Thus  in  this  case 
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the  laminate  is  assumed  to  be  stress  free  at  room  temperature. 
This  difference  in  internal  states  at  room  temperature  is 
reflected  by  the  difference  in  measured  thermal  expansion  co¬ 
efficients  for  these  two  cases.  Results  of  analytical  calcu¬ 
lations  based  upon  the  above  assumptions  are  shown  in  the 
lower  section  of  table  5.  The  same  trend  is  calculated  as 
was  observed.  For  example,  for  a  value  of  of  69  MPa,  the 
calculated  thermal  expansion  coefficient  is  1,17;  while  when 
the  delay  is  assumed  to  result  in  a  stress  free  temperature 
equal  to  room  temperature,  the  calculated  expansion  coeffi¬ 
cient  for  the  subsequent  cooling  cycle  is  equal  to  2.27,  As 
can  be  seen  in  tables  3c  and  3d  the  respective  experimental 
numbers  are  1.17  and  4.74.  The  expansion  coefficients  in  the 
RT  to  -54°C  range  increase  when  there  is  a  long  delay  at  room 
temperature.  However,  the  measured  increase  is  even  larger 
than  the  calculated  one. 

An  additional  result  of  interest  is  the  significant  dif¬ 
ference  in  thermal  expansion  coefficients  for  the  cooling  and 
heating  parts  of  the  cycle  RT  -54°C  ^  RT  as  shown  in  table 
3d.  This  can  be  explained  as  follows:  it  is  seen  that  the 

t 

internal  stress  0^2  attained  at  -54°C  is  a  high  tension 
(table  5,  RT  ^  -54°C,  no  delay  at  RT) .  Thus  the  thermal  ex¬ 
pansion  coefficient  in  this  range  is  significantly  affected 
by  material  nonlinearity.  Upon  reheating  from  -54 °C  to  RT,  the 

I 

transverse  stress  a unloads  and  it  may  be  considered  that 
the  internal  strains  diminish  elastically .  Therefore,  the 
thermal  expansion  coefficient  in  this  heating  range  is  one 
for  an  elastic  laminate  with  final  elastic  properties  those 
of  room  temperature.  For  this  temperature  range  that  elastic 
expansion  coefficient  is  1.51  as  compared  to  the  calculated 
nonlinear  value  for  this  temperature  range  of  1.17  (utilizing 
=  69  MPa) .  The  experimental  result  (see  table  3d)  for 
heating  over  this  temperature  range  is  larger  than  for  cool¬ 
ing.  This  agrees  qualitatively  with  the  calculated  result, 
but  the  quantitative  difference  is  much  larger  than  calcu¬ 
lated.  This  suggests  the  significance  of  nonlinearities  and 


time  dependence  of  the  material. 

Influence  of  Applied  Loads 

The  next  set  of  tests  is  concerned  with  thermal  expan¬ 
sion  of  a  ±45°  laminate  when  subjected  to  uniaxial  compressive 
loads.  Table  6  lists  the  measured  values  of  the  strains  in 
the  load  direction  and  perpendicular  to  the  load  for  various 
temperatures.  In  the  first  part  of  the  table  the  applied 

O 

stress  was  a..  =  35.4  MPa,  while  in  the  second  the  applied 

oil  I 

stress  was  =  70.8  MPa.  The  shear  stress 

fibers  in  ±45°  direction  then  follows  from  equation  (2.3.5) 

I 

and  the  corresponding  shear  strain  follows  from  equation 

(2.3.12)  as: 

I  O  O 

^^12  "  ^11  "  ^22 

and  is  listed  in  table  6. 

Values  of  the  shear  modulus  at  various  temperatures  have 
been  extracted  from  the  data  using  the  ratios  between  the  con¬ 
stant  shear  stress  and  the  measured  shear  strain  at  that  tem- 

o 

perature .  The  loading  =  -35.4  MPa  has  been  used  since 

it  can  be  assumed  that  for  such  low  loading  nonlinearities  in 
shear  stress-strain  relations  are  not  significant.  The  shear 

o 

moduli  for  the  loading  ~  "70.8  MPa  are  smaller  than  for 

the  lower  loading.  This  can  be  explained  by  lamina  nonlinear¬ 
ity  which  becomes  more  significant  at  higher  loading. 

Plots  of  lamina  strains  versus  temperature  for  the  var¬ 
ious  loadings  are  shown  in  figures  5-8.  The  behavior  is  non¬ 
linear  and  there  is  a  large  reduction  in  material  properties 
at  temperatures  beyond  149°C.  This  accounts  for  the  increases 
in  strain  in  this  temperature  range. 

Experimental  results  for  the  thermal  expansion  coeffi¬ 
cients  for  the  loaded  laminate  are  shown  in  table  4.  It  is 
seen  that  these  expansion  coefficients  are  quite  different 
from  the  free  thermal  expansion  coefficients  of  the  laminate, 
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table  3.  Furthermore,  the  expansion  coefficients  transverse 
to  the  load  direction  are  quite  different  from  those  in  the 
load  direction. 

Tables  7-10  summarize  computation  of  thermal  expansion 
coefficients  for  laminates  loaded  in  compression  on  the  basis 
of  the  theory  of  paragraph  2.3.  The  results  are  given  for 
three  loadings:  ^ ~  -35.4;  -70.8;  -106.  MPa;  and  for 

various  combinations  of  "yield  stresses"  for  transverse 
stress,  a  ,  and  for  shear  stress,  t  .  For  the  elastic  case 

y  y  *  * 

(cr  =  T  =  co)  the  expansion  coefficients  and  a2  the 

same  and  equal  to  the  free  thermal  expansion  coefficients  of 
the  laminate  (see  table  5) .  (Small  differences  are  due  to  dif¬ 
ference  in  properties  at  various  temperatures  considered.) 

The  calculated  results  show  that  nonlinearities  have  a 

significant  influence  on  the  thermal  expansion  coefficients. 

*  .  * 

Decrease  of  a  significantly  decreases  and  increases  a^. 

Decrease  of  x  has  a  similar  but  much  less  significant  effect, 

y 

As  an  example,  the  results  for  =  69,  and  x^  =  138  MPa 
illustrate  the  effects  of  nonlinearity.  In  the  temperature 
range  RT  to  93. 3*^0:  =  1.86  and  CL2  ~  3.78  (in  units  of 

y  cm/cm/°C) ;  while  for  an  elastic  laminate  =  3.01. 

Thus,  the  theory  predicts  that  the  nonlinear  loaded  laminate 
is  anisotropic  in  thermal  expansion  and  its  thermal  expansion 
coefficients  are  not  equal  to  the  free  thermal  expansion  co¬ 
efficients.  The  experimental  results  show  the  same  effect. 

Discussion 

While  the  experimental  results  agree  qualitatively  with 
the  analytical  predictions,  the  actual  numerical  agreement 
cannot  be  considered  satisfactory.  There  are  several  possible 
reasons  for  the  discrepancies. 

1 .  Measured  values  of  thermal  expansion  coefficients  of 
the  unidirectional  material  have  shown  considerable 
scatter,  in  particular  for  the  longitudinal  thermal 
expansion  coefficient  a^.  Lamina  thermal  expansion 
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coefficients  used  in  laminate  analysis  were  averages 
of  measured  thermal  expansion  coefficients  for 
several  different  specimens. 

2.  The  samples  supplied  for  property  measurements,  both 
unidirectional  and  ±45°  laminates,  were  not  of  high 
quality.  Examination  under  microscope  revealed  many 
defects  and  flaws. 

3.  The  analysis  was  based  on  the  assumption  that  the 
unidirectional  material  is  nonlinear  in  response  to 
transverse  and  shear  stresses  for  loading  and  un¬ 
loading.  A  more  realistic  approach  may  be  to  con¬ 
sider  nonlinearity  in  loading  and  conventional 
elastic  behavior  in  unloading.  Indeed,  one  striking 
effect  in  thermal  expansion  behavior  (the  significant 
difference  in  expansion  coefficient  in  the  two 
branches  of  the  thermal  cycle  RT  -54°C  RT)  could 
be  explained  by  elastic  unloading. 

4.  Thermal  expansion  analysis  presupposes  a  stress  free 
state  at  some  known  reference  temperature.  In  a  non¬ 
linear  laminate  the  values  of  the  thermal  strains 
depend  on  the  value  of  the  stress  free  temperature. 
However,  the  stress  free  temperature  is  only  imper¬ 
fectly  known  since  it  depends  on  previous  history  of 
the  laminate.  The  importance  of  this  phenomenon  has 
been  demonstrated  here  by  the  differences  in  free 
thermal  expansion  coefficients  due  to  delay  times  at 
room  temperature  when  a  laminate  is  cooled  from 

260® to  -54°C.  It  may  be  assumed  that  a  laminate  is 
stress  free  at  its  curing  temperature.  When  it  is 
cooled  to  room  temperature  significant  thermal 
stresses  develop.  These  stresses  relax,  thus  dimin¬ 
ish  with  time  at  room  temperature.  To  assess  this 
phenomenon  quantitatively  a  detailed  analysis  of  the 
time-dependent  relaxation  of  stresses  in  a  laminate 
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is  needed.  This  requires  measurement  of  temperature 
dependent  relaxation  moduli  and  creep  compliances  of 
a  lamina  and  their  incorporation  into  laminate  analy¬ 
sis.  A  more  fundamental  approach  would  be  to  measure 
temperature  dependent  relaxation  moduli  and  creep 
compliances  for  the  matrix  and  then  compute  such 
quantities  by  micromechanics  analysis  for  the  uni¬ 
directional  material  (see  [3]). 


CONCLUSIONS 


The  effects  demonstrated  in  this  study  have  serious  implica¬ 
tions  for  analysis  of  laminates  which  behave  nonlinearly  in  the 
presence  of  significant  temperature  changes.  One  of  the  most 
basic  engineering  properties  of  a  laminate  is  its  thermal  expan¬ 
sion  coefficient  which  determines  the  thermal  strain  due  to  unit 
temperature  rise.  When  a  thermoelastic  laminate  is  heated  and 
loaded  the  strains  due  to  load  and  heating  are  simply  superposed. 
This  procedure  is  still  valid  when  the  linear  thermoelastic  pro¬ 
perties  of  the  laminae  are  temperature  dependent.  It  has  been 
demonstrated  in  this  study,  both  analytically  and  experimentally, 
that  this  superposition  becomes  invalid  when  the  mechanical  pro 
perties  of  the  laminae  are  nonlinear  -  that  is:  the  stiffnesses 
are  dependent  on  the  state  of  stress.  In  this  case  the  free 
thermal  expansion  coefficient  is  not  useful  to  compute  the  ther¬ 
mal  strains  of  a  laminate  due  to  temperature  change  in  the  pres¬ 
ence  of  applied  loads.  The  thermomechanical  strains  and  the  in¬ 
ternal  stresses  must  be  computed  by  thermomechanical  nonlinear 
laminate  analysis. 

In  summary: 

1.  Material  nonlinearity  has  a  significant  effect  on  the 
thermal  expansion  coefficients  of  laminates. 

2.  The  thermal  expansion  coefficients  for  an  unloaded  and  a 
loaded  laminate  are  significantly  different  when  the 
laminate  deforms  in  a  nonlinear  manner.  For  such  a 
loaded  laminate  the  thermal  expansion  coefficients  are 
different  in  direction  of  load  and  transverse  to  load. 

3.  Experimental  data  obtained  for  a  symmetric  0/90  laminate 
indicate  that  time  dependence  of  the  material  and  the 
associated  stress  relaxation  significantly  affect  the 
values  of  free  thermal  expansion  coefficients. 
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of  Thermal  Expansion  Coefficients  with  Temperature 


Table  3.  Effect  of  Temperature  and  Time  Delay  on  Coefficients 
of  Free  Thermal  Expansion  for  ±45°  Laminates, 

(ye/^C) 


(a)  Thermal  Expansion  without  Time  Delay 


Sample 

Ascending 

10°  -  149°C 

Descending 

149°  -  10°C 

D5 

3.87 

5.17 

C2 

5.09 

4.75 

C3 

4.27 

5.42 

Mean 

4.41 

5.11 

(b)  Thermal  Expansion  with  Time  Delay 


Sample 

260°C  RT 

Delay  Time 

RT  -*■  260°C 

D4 

3.13 

15  Days 

3.29 

C4 

3.24 

15  Days 

5.36 

D6 

1.80 

15  Days 

2.52 

Mean 

2.74 

15  Days 

3.73 

Table  3.  Continued 


(c)  Thermal  Expansion  with  Time  Delay 


Sample 

260°C  RT 

Delay  Time 

C5 

2.12 

15  Days 

3.13 

D7 

1.4  8 

15  Days 

6.19 

Cl 

4.32 

15  Days 

4.91 

Mean 

1.47 

15  Days 

4.74 

(d)  Thermal  Expansion  without  Time  Delay 


Sample 

260‘=’C  RT 

Delay  Time 

RT  -54°C 

-54‘^C  RT 

C5 

2.12 

None 

.52 

5.06 

D7 

1.48 

None 

.  70 

6.00 

Cl 

4.32 

None 

2.29 

7.16 

Mean 

2.65 

None 

1.17 

6.34 
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of  Thermal  Expansion 
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Table  5,  Calculated  Thermal  Expansion  Coefficients  of 

G/PI  ±45®  Laminate  (No  externally  applied  loads) 
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Table  6.  Shear  Strains  on  Loaded  ±45“  Laminate 


♦ 

0 

^22 

0 

^11 

O 

n  .  =  -35.37  MPa 

11 

-51 

2.025 

-3.185 

2.605 

-17.72 

-18 

2.150 

-3.050 

2.600 

-17.72 

10 

2.285 

-2.900 

2.539 

-17.72 

22.2  (RT) 

2.395 

-2.800 

2.598 

-17.72 

38 

2.535 

-2.675 

2.605 

-17.72 

66 

2.685 

-2.525 

2.605 

-17.72 

93 

3.150 

-2.475 

2.813 

-17.72 

149 

3.510 

-2.400 

2.955 

-17.72 

204 

3.685 

-2.700 

3.193 

-17.72 

Oil  =  70.74  MPa 


-51 

-18 

10 

22.2  (RT) 


4.550 

-6.275 

5.413 

35.37 

4.750 

-6.100 

5.425 

35.57 

4.975 

-5.700 

5.338 

35.37 

5.120 

-5.590 

5.355 

35.37 

5.300 

-5.450 

5.375 

35.37 

5.800 

-5.450 

5.625 

35.37 

6.300 

-5.700 

6.000 

35.37 

7.050 

-6.550 

6.800 

35.37 

8.200 

-8.300 

8-.  250 

35.37 

<|)  (Temperature)  -  “C 
Stress  ~  MPa 

Strain  “  10 


-  10"^  SS 

cm 


Table  7.  Effect  of  Temperature  and  Load  on 

Calculated  Thermal  Expansion  Coefficients 
and  Stresses  of  ±45'*  G/PI  Laminate 
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Table  9.  Thermal  Expansion  Coefficients  and  Internal  Stresses 
of  ±45“  Loaded  Laminate 
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A  =  Active  gage  mounted  on  polyimide 
C  =  Compensation  gage  mounted  on  quartz 


Figure  4.  Half  Bridge  Circuit  Used  in  Temperature  Compensated 
Tests 
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Figure  7.  LongitudJLn»J.  Strain  versus  Temperature  Load 
=  -2946  M  (-300  kg)  =  -70.74  MPa 
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